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The spatial curvature (K or Ll K ) is one of the most fundamental parameters of an isotropic and homogeneous 
universe and has a close link to the physics of the early Universe. Combining the radial and angular diameter 
distances measured via the baryon acoustic oscillation (BAO) experiments allows us to unambiguously constrain 
the curvature. The method is primarily based on the metric theory, but is less sensitive to the theory of structure 
formation other than the existence of the BAO scale and is free of any model of dark energy. In this paper, we 
estimate a best achievable accuracy of constraining the curvature with the BAO experiments. We show that an 
all-sky, cosmic-variance-limited galaxy survey covering the Universe up to z i 4 enables a precise determination 
of the curvature to an accuracy of cr(fl K ) ^ 10 -3 . When we assume a model of dark energy - either the 
cosmological constraint or the (wo,w a ) model - it can achieve a precision of cr(Q. K ) ^ a few x 1CT 4 . These 
forecasts require a high sampling density of galaxies, and are degraded by up to a factor of a few for a survey 
with a finite number density of ~ KT 3 (/i/Mpc) 3 . 

PACS numbers: 98.80.-k,95.36.+x,95.75.-z,98.65.Dx 

I. INTRODUCTION 

The curvature of the Universe (hereafter denoted as K or 
Q f () is one of the most fundamental quantities in an isotropic 
and homogeneous universe in the context of general relativ¬ 
ity (GR) (T). The curvature also has a close connection to 
the physics of the early Universe. An inflationary universe 
scenario predicts that the “apparent” curvature, which we can 
infer from an observable universe, should appear to be close 
to a flat geometry (Qx - 0), even if the exact value is nonzero 
01 HI- If the Universe arose from the decay of a false vacuum 
via quantum tunneling, it leads to an open geometry (K < 0 
or L1 k > 0) BE). In particular, if the Universe began with 
“large-field inflation” @ - which predicts that the primordial 
gravitational wave is as large as can be observed in the cos¬ 
mic microwave background (CMB) anisotropies - or if the 
universe began with “just enough inflation” in a landscape or 
multiverse picture, the curvature can be large enough to be 
measurable (Q K ~ 10 4 —10 2 ) liT0UT4l . Further, an addition 
of the super-curvature perturbation in an open-inflation sce¬ 
nario might resolve the large-scale CMB anomalies D3H3- 
On the other hand, if a closed curvature (K > 0 or Qx < 0) is 
found, it gives rise to a challenge for the inflationary scenario: 
the Universe needs to emerge from a specific initial condition 
of closed curvature 13 ED. Thus it is important to constrain 
the curvature from cosmological observations in order to ob¬ 
tain a clue to the physics of the early universe [also see Ref. 
m for a thorough review of various cosmological probes]. 

Observations of the CMB have led to the precise measure¬ 
ment of the angular diameter distance to the last scattering sur¬ 
face, preferring a flat geometry as expected in an inflationary 
scenario [e.g.,[20l[2T|- However, the constraint rests on the use 
of the standard cosmological model such as the cosmological 
constant, cold-dark-matter dominated model (ACDM model), 
based on GR and the nearly adiabatic initial conditions. If 


the assumptions are relaxed - for instance, if a generalized 
model of dark energy is employed - the CMB constraint on 
the curvature is largely degraded (22). The baryonic acoustic 
oscillation (BAO) provides us with an alternative, powerful 
geometrical probe, allowing one to constrain the cosmologi¬ 
cal distances via measurements of the galaxy clustering pat¬ 
tern in redshift and angular directions, respectively 1 12314261 . 
The BAO experiments are shown to be robust against various 
astrophysical systematic effects such as the galaxy bias un¬ 
certainty |[23[28ll . The current state-of-the-art measurements 
were done using the Sloan Digital Sky Survey III Baryon Os¬ 
cillation Spectroscopic Survey (291 - 1311 . achieving a percent 
precision of the distance measurement at z — 0.57. Various 
wide-area spectroscopic galaxy surveys are planned that aim 
to achieve precise BAO measurements up to higher redshifts: 
the Subaru Prime Focus Spectrograph (PFS) (32), the Dark 
Energy Spectrograph Instrument (DESI) (351 . the ESA Eu¬ 
clid satellite mission (57), and the NASA WFIRST mission 

ea. 

The BAO method is unique in that it can constrain the radial 
(more exactly the Hubble expansion rate) distance as well as 
the angular diameter distance at the redshift of the galaxy sur¬ 
vey, while other geometrical probes such as supernovae and 
gravitational lensing can probe the luminosity or angular dis¬ 
tances (and not the radial distance). The relation between the 
radial and angular diameter distances is purely geometrical 
and specified by the curvature; if the Universe has a nonzero 
curvature, the two distances differ. The relation holds for any 
theory of gravity or dark energy, and rests on the metric theory 
of a homogeneous and isotropic space, which is a maximally 
symmetric spacetime described by the Friedmann-Robertson- 
Walker (FRW) metric m. Hence the purpose of this paper is 
to estimate the fundamental accuracy of estimating the curva¬ 
ture parameter with the BAO experiments [see Refs. [35} [36} 
for the previous works based on the similar motivation]. To 
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do this, we will assume a cosmic-variance-limited galaxy sur¬ 
vey, namely full-sky coverage and a sufficiently high number 
of density of sampled galaxies. 

The structure of this paper is as follows. In Sec.[EI]we dis¬ 
cuss the methodology based on the Fisher matrix information 
formalism, and briefly review the BAO method. In Sec. [HI] 
we show the main results. Section m is devoted to discus¬ 
sion. Unless stated otherwise, we will adopt as fiducial model 
a flat ACDM cosmology with £2 m o = 0.27, Qa = 0.73 and the 
Hubble parameter h = 0.71. 


II. GEOMETRICAL ESTIMATION OF THE CURVATURE 
WITH BAO DISTANCES 


A. Cosmological distances 


We assume that the Universe is statistically isotopic and 
homogeneous. The spacetime structure of such an universe 
is described by the FRW metric m. With the metric theory, 
solving the light propagation in an expanding universe yields 
a relation between cosmological distances and redshift. The 
comoving radial distance is given in terms of the integral of 
the Hubble expansion rate: 


O c (z) 


f 


dz' 


o Hiz'Y 


(1) 


where the Hubble expansion rate is given by the time deriva¬ 
tive of the scale factor as H = din a/dt, and z is the redshift, 
given as 1 + z = 1/a where we employed the convention 
a (to) = 1 today. The metric theory also gives a geometrical 
relation between the radial and angular diameter distances: 


D a (z) = 


_ _ sinh Y^KDc(z) (K < 0), 

\Tk 

D c {z) (K = 0), 

-2= sin YkD c (z ) (K > 0) 

Vk 


( 2 ) 


The different equations are for different geometries of the Uni¬ 
verse; open, flat and close geometries, respectively. The cur¬ 
vature K is in units of 1/distance 2 . The CMB experiments 
[e.g.. Ref. |20j| give stringent constraints on the curvature, im¬ 
plying = \K\f if- <s 1, where Q. K is the curvature param¬ 
eter and 1/Ho is the Hubble radius today. Note that through¬ 
out this paper we use the comoving angular diameter distance, 
rather than the physical distance, and the two are related as 
/J phys (z) = Da (z) l (1 -i-s). For redshifts relevant for galaxy sur¬ 
veys, where y/\K\Dc(z) « 1 holds, Eq. j2j can be expanded 
around \l\K\D ( (z) = 0, yielding the approximation 


D a (z ) - D c (z) 


1 - l -KD c (zf 
o 


(3) 


Both D A (z) and Dc(z) are observables of the BAO experi¬ 
ments for each redshift slice. In the following we will use 
the above equation to estimate the accuracy of estimating the 
curvature from BAO information without making many as¬ 
sumption about the theory of structure formation (other than 


the existence of the BAO scale in the distribution of galaxies) 
or any model of dark energy. 


B. Estimator of the curvature and its covariance 


Suppose we have the BAO distance measurements in N s 
redshift bins, without any gap, over a range of redshift from 
today up to a maximum redshift z max : D A (zi ) and DhIz/) for 
i = 1,2,...,N S (zi = 0 and zn, = Zmax)- Here D H (z) is the 
comoving Hubble distance: Dh(z ) = 1 /I Hz.). 

The radial distance at the ith redshift bin n can be estimated 
by combining the measured Hubble distances over a redshift 
range z = [0, z,] (see Eq. |T]i: 

D c ,i - 2 D H jAzj, (4) 

Zj<Zi 

where we have introduced the notations, D ci = Dcizj and 
D Hi = I) H (zi). In practice one might want to use a more so¬ 
phisticated method to estimate £>c(z) to avoid inaccuracy due 
to the discrete summation, e.g. by assuming that H(z) is mod¬ 
eled by a polynomial function of redshift and then estimating 
the coefficients from fitting of the function to the measured 
H(z). Here we assume a discrete summation for simplicity. 

From Eq. ([3]), an estimator of the curvature parameter from 
measurements of the radial and angular diameter distances at 
the ith redshift bin can be given as 


K, = 6 


D A ,i - £>c,i 

Dh 


(5) 


Combining the measurements at different redshift bins, we 
can define the chi-square (x 2 ) to estimate the curvature pa¬ 
rameter as 


X 


= K)\C kk \' ( Kj-K ), 


Uj= 1 


( 6 ) 


where K is the underlying true curvature, treated as a free 
model parameter in the x~ fitting. [C KK ] 1 denotes the in¬ 
verse of the covariance matrix C KK defined as 


r<KK 
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Here the covariance matrices such as C DcDc describe statis¬ 
tical uncertainties of the BAO observables - the Hubble and 
angular-diameter distances - given BAO measurements of a 
galaxy survey: 

C? cDc = Cov[D c „ Dcd = XX C'Yf H &z m Az,u 

m<i n<j 

C° cD * = Cov[D Ci , bed = Z C lf*Az m , (8) 

m<i 
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and so on. Around a flat-geometry universe, where Da — D c , more quantitative, the redshift-space galaxy power spectrum 
the covariance matrix is approximated as at redshift z is given in the linear regime as 


D A,i D AJ KK 

36 ij 


* Cf. 


- [c° aDc + C y . aD ' I + C® cDc (9) 


The statistical uncertainty of the curvature estimation is 
given from the second derivatives of x 2 (K) (Eq. [6} with re¬ 
spect to the true curvature: 


cr 2 (K) = 



( 10 ) 


where \ C KK ] ! is the inverse of Eq. ([7]». We will use the above 
equation to estimate the accuracy of the curvature estimation 
for a given galaxy survey. 


C. BAO 


In this subsection, assuming a hypothetical galaxy survey, 
we derive the Fisher information matrix of the Hubble and 
angular diameter distances from the BAO measurements. 

The two-point correlation function or the Fourier- 
transformed counterpart - the power spectrum - is measured 
as a function of the separation lengths between paired galax¬ 
ies. In this procedure, the position of each galaxy needs to 
be inferred from the measured redshift and angular position. 
Then the separation lengths perpendicular and parallel to the 
line-of-sight direction from the measured quantities are given 
as r i oc A 9 and ry oc Az, where A 6 and Az are the differences 
between the angular positions and the redshifts of the paired 
galaxies. For this conversion, we need to assume a reference 
cosmological model to relate the observables ( Ad , Az) to the 
quantities (rj_,ry). Thus, the wave numbers are given as 


^.L.ref — 


D a (z ) , 

D A sef(z) 


^||, ref — 


Dh(z) 

Dh.kKz) 


k\\. 


( 11 ) 


The quantities with the subscript “ref” are the quantities es¬ 
timated from the observables assuming a “reference” cos¬ 
mological model, and the quantities without the subscript 
are the underlying true values. Since the reference cosmo¬ 
logical model assumed generally differs from the underlying 
true cosmology, it causes an apparent distortion in the two- 
dimensional pattern of galaxy clustering. In principle, this 
could be measured using only the isotropy of clustering statis¬ 
tics - the so-called Alcock-Paczynski (AP) test Ii37l - but a 
more robust measurement of both Da(z) and Dh(z) can be ob¬ 
tained by searching for the “common” BAO scales in the pat¬ 
tern of galaxy clustering, as the standard ruler, in combination 
with the CMB constraints Il23lf24l[38ll . 

We will use the currently standard ACDM model as a guid¬ 
ance for the parameter dependence of our constraints and as 
an effective realistic description of the galaxy clustering. Nev¬ 
ertheless, we would again like to emphasize that the method¬ 
ology proposed in this paper relies only on an FRW metric 
theory and the existence of the BAO standard ruler. To be 


Pg,s(k_ L.ref? ^||,ref j z) 


Dam{z) Oftreffe) 
D a (z) 2 D h (z) 


k 2 1 2 

1 +p{z) i 


b]P L m {t,z) + Psn,{ 12) 


where b g is the linear bias parameter, /j is the lin¬ 
ear redshift-space distortion parameter, defined as [3 = 

(1 /b g ) d\nD/d\na\ z l39l . D is the linear growth rate, P Y (k) is 
the linear mass power spectrum, and P sn is a constant number 
and is a parameter to model the residual shot noise. Through¬ 
out this paper we assume b(zi) = 1 as the fiducial value in each 
redshift bin for simplicity. This is a conservative assumption, 
because most galaxies at higher redshift are very likely to be 
biased tracers with b > 1. When deriving the BAO geomet¬ 
rical constraints we will marginalize over the effect of galaxy 
bias uncertainty. 

To make the parameter forecast, we employ the method de¬ 
veloped in Ref. l28l . In this method, we include the smearing 
effect of the BAO features due to the bulk flow of galaxies in 
large-scale structure t4QH42|. For the BAO survey of multiple 
redshift bins, the Fisher information matrix of model parame¬ 
ters can be computed as 


pgalaxy 
r afj 


X 


f 2 n)cdk 

Jw w 

d In Pgjk, p\zi ) d In P g , s (k, p\ zd 
dp a dpp 



xV eff (k; Zi ) exp [-k 2 !, 2 ± - k 2 p 2 (^ - 5*)], 


(13) 


where p is the cosine between the wave vector and the line-of- 
sight direction, p = fcy/k; 2- is the sum over different redshift 
bins; dP gtS /dp a is the partial derivative of the galaxy power 
spectrum (Eq.[l2| with respect to the crth parameter around the 
fiducial cosmological model; and the effective survey volume 
y eff and the Lagrangian displacement fields Ey and £ to model 
the smearing effect are given as 


Ves(k, p; zd 

X±(z) 

2||(z) 


n g (Zi)P g , s (k, P', zd 
n g (Zi)PgJk,p\Zi) + 1 


Fsurvey(^f)? (14) 


c rec £»(z)£o, (15) 

c iec D(z)(l +/ 9 )£o- (16) 


Here V S urveyfo) is the comoving volume of the redshift slice 
centered at z,; the present-day Lagrangian displacement field 
is £o = 11 /z^'Mpc for cr 8 = 0.8 l27ll ; D(z) is the growth rate 
normalized as D(z = 0) = 1; f g — d In D/d In a. The param¬ 
eter c rec is a parameter to model the reconstruction method of 
the BAO peaks (see below). In Eq. ( fT3j ), we take the expo¬ 
nential factor of the smearing effect outside of the derivatives 
of P g s . This is equivalent to marginalizing over uncertainties 
in £y and £ ± . The growth rate in £y or £ ± takes into account 
the smaller smearing effect at higher redshift due to the re¬ 
duced evolution of large-scale structure. For the parameters, 
we include the cosmological parameters, the distances in each 
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redshift slice, and the nuisance parameters: 

Pa — '4 v, n s , (y s . Q nl (,/r, fT|/p, k) A (Zi). 

O w (z,), bg(zd,P(zd, /’safe)}, (17) 

where A v , /z v and a s are parameters of the primordial power 
spectrum, A s is the amplitude of the primordial curvature per¬ 
turbation, and n s and a s are the spectral tilt and the running 
spectral index. The set of cosmological parameters deter¬ 
mines the shape of the linear power spectrum. For the k in¬ 
tegration, we set k m [„ = 1CF 4 h /Mpc and k max = 0.5 /i/Mpc 
for all the redshift slices, but the exponential factor in Eq. ( | 1 3[ > 
suppresses the information from the nonlinear scales. The 
Fisher parameter forecasts depend on the fiducial cosmolog¬ 
ical model for which we assumed that the model is consis¬ 
tent with the WMAP 7-year data l20l . For a galaxy survey of 
N s redshift bins, the number of model parameters is in total 
6 + 5 x N s . 

Furthermore, we assume the BAO reconstruction method 
in Ref. 127). Since the peculiar velocity field of galaxies in 
large-scale structure can be inferred from the measured galaxy 
distribution, the inferred velocity field allows for pulling back 
each galaxy to its position at an earlier epoch and then recon¬ 
structing the galaxy distribution more in the linear regime. As 
a result, one can correct to some extent the smearing effect 
in Eq. 0 and sharpen the BAO peaks in the galaxy power 
spectrum. Padmanabhan et al. Il43l implemented this method 
to the real data, SDSS DR7 LRG catalog, and showed that 
the reconstruction method can improve the distance error by 
a factor of 2. The improvement was equivalent to reducing 
the nonlinear smoothing scale from 8.1 to X nl = 4.4 h ' .VIpc, 
about a factor of 2 reduction in the displacement field. To 
implement this reconstruction method requires a sufficiently 
high number density of the sampled galaxies in order to reli¬ 
ably infer the peculiar velocity field from the measured galaxy 
distribution [see Ref. [30] for the latest result]. In the Fisher 
matrix calculation, we used c rec =0.1 for an implementation 
of the reconstruction method l58l . 

The BAO reconstruction has so far been successful in sub¬ 
stantially improving the accuracies of both the angular diam¬ 
eter distance and radial Hubble distance estimations [see Ref. 
[441 for a recent review]. The systematic errors can arise from 
an imperfect treatment of nonlinear effects such as nonlinear¬ 
ities in the density field and the nonlinear motions of galax¬ 
ies inside massive halos. A further refinement of the recon¬ 
struction method will be required, by using mock catalogs 
of galaxy surveys, in order to achieve the ultimate precision 
[e.g.,04] for such an attempt]. A striking advantage of using 
a higher redshift galaxy survey is that the nonlinear effects are 
weaker at higher redshifts where large-scale structure evolves 
less, and it will therefore allow a more robust reconstruction 
to obtain unbiased measurements of the BAO distances. 

In the following forecast, we assume the BAO experiments 
combined with the CMB constraints expected from the Planck 
satellite: 

F = F cmb + F lgalaxy , (18) 

where Fcmb is the Fisher matrix for the CMB measurements. 
We employ the method in Ref. Il32ll to compute the CMB 


Fisher matrix, where we assume the standard ACDM model 
for the physics prior to recombination that determines the 
sound horizon scale or the BAO scale. 

To estimate the accuracy of the distance estimation with 
the BAO experiments, we first invert the Fisher matrix F and 
then use the submatrix including elements of the Hubble and 
angular diameter distances: 

invert , use submatrix , 

F => [FT 1 => [F]~mn (19) 

where the indices M,N denote the elements including the 
Hubble and angular diameter distances. The submatrix Fmn 
gives the error covariance matrices of distances, D A (z ,■) and 
Dn(Zi), after marginalizing over other parameters (6 cosmo¬ 
logical parameters and nuisance parameters). Hence the di¬ 
mension of Fmn is (2 N s ) x (IN,). The covariance matrices of 
D A {zd and D H {Zi) and the cross-covariance matrix are 

C-/ Da = [f-'] D a(Zi)Da(Zj) , 

C D A D H _ r ri 1 "I 

ij = \D A (Zi)D H ( Z j)> 

C°" D " = [F- 1 ] Dh(z , )Dh(Zj) . (20) 

We use these covariances in Eq. 0 to estimate the accuracy 
of the curvature estimation. 


HI. RESULTS 
A. Survey parameters 

Here we assume an ideal survey to estimate the fundamen¬ 
tal limit of the curvature estimation via the BAO measure¬ 
ments. That is, we assume a galaxy survey with full-sky 
coverage and a sufficiently high number density of sampled 
galaxies in each redshift bin up to a given maximum redshift 
z max , without any gap. We will study how the curvature deter¬ 
mination is degraded when considering a survey with a finite 
number density of galaxies. 


B. Forecasts 

Figure [I] shows the forecasts of the distance measurements 
with a cosmic-variance-limited survey up to z max = 7. The 
shaded region around each redshift bin with width Az = 0.2 
denotes the +lcr uncertainty in the angular diameter distance 
including marginalization (Eq.|20)i. The accuracy of the Hub¬ 
ble distance is only slightly worse than the error in D,\ (at the 
level of 10%). As described in Sec. |IIC| we include both the 
BAO feature and the broadband feature of the galaxy power 
spectrum (the AP test) for the distance forecasts. We con¬ 
firmed that the distance accuracies are mainly from the BAO 
features. We should also stress that, once a sufficient number 
of redshift slices, up to z ~ 1 for our setting, are included, the 
BAO information of the galaxy survey allows a more accu¬ 
rate determination of the sound horizon scale than the Planck 
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FIG. 1 : The shaded region in each redshift bin, with width Az = 0.2, denotes the fractional error of the angular diameter distance determination 
D A (zi), expected for a cosmic-variance-limited galaxy survey with full-sky coverage (/ S k y = 1) and a sufficiently high number density of 
galaxies in each redshift over the redshift range from z = 0 to z max = 7. The errors include marginalization over other parameters (see text 
for details). For comparison, the solid curves show the fractional difference in the radial or angular diameter distances, D c (z) or D A (z), for an 
open-geometry universe with fl K = 0.001 compared to the fiducial flat universe, with the dark matter density and the cosmological constant 
being fixed. The curvature parameter alters D c (z) and D A in different ways and the ratio, D A /D C , displays a characteristic redshift dependence. 
To be more precise, the fractional difference is defined, e.g. as D A (z,Cl K = 0.001)/Z)A,fid(z; Qk = 0) — 1. The dashed and dot-dashed curves 
show the fractional difference when changing n m0 and the dark energy equation of state w 0 , by A£2 m0 = 0.0001 and Aiuo = 0.01, within a flat 
universe. In this case, the changes keep D A (z) = D c (z). 


prior, which explains why the distance accuracies get satu¬ 
rated at z ~ 1 ■ 

Using the forecasts of the distance measurements, we show 
in Fig. [2] the expected accuracy of the curvature estimation, 
ct(Q. k ). Here we consider a cosmic-variance-limited galaxy 
survey over a range of redshift z = [0, z max ] and show the 
forecasts as a function of the maximum redshift z max . Note 
that the error cr(Q. K ) is free of the uncertainty in the Hubble 
parameter h, because the uncertainty in h is absorbed when 
estimating the curvature parameter from a combination of the 
radial and angular distances. The top curve shows the fun¬ 
damental limit, estimated from Eq. (|7j, in the sense that the 
constraint is purely based on the BAO distance measurements 
and is less sensitive to the theory of structure formation (other 
than the BAO scale) such as the galaxy bias and is free of 
any model of dark energy. The accuracy can be better than 
0.1%, if a cosmic-variance-limited survey up to z max i 5 is 
available. Since the curvature has a relatively larger impact 
on the expansion history up to higher redshifts than the cos¬ 
mological constant, adding the BAO distance measurements 
at higher redshift keeps improving the curvature constraint. 

If we employ Einstein gravity (GR), which relates the struc¬ 


ture of spacetime to the energy-matter content of the Universe 
through the Einstein equations, the curvature constraint can 
be improved, although the constraint is now model-dependent. 
At redshifts relevant for galaxy surveys, the Hubble expansion 
rate is given as 

H(z) 2 = Hq fo m0 ( 1 + z) 3 + a K (1 + zf + £2 de (z)] • (21) 

Then the Hubble and angular-diameter distances are specified 
by the density parameters (see Eqs. [T] and [2j. However, the 
nature of dark energy is another mystery of the universe. For 
comparison, we here employ a phenomenological model of 
dark energy that is parametrized in terms of its equation-of- 
state parameters as 

wde(z) - w 0 + w a (l - a), (22) 

where wq and w a are parameters and constant in time. In 
this case, the redshift evolution of dark energy is given as 
PDe(z) a -3(i+<Do+!0o) e -3<o,,(i-a) w jth 1 + z — 1 /a. We should 
emphasize that this model acts as a “strong” prior, restricting 
the degrees of freedom of the dark energy model to a spe¬ 
cific model. The redshift dependence of the above dark en¬ 
ergy model, around the fiducial cosmological constant model. 
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FIG. 2: An expected accuracy of the curvature parameter determination with the BAO measurements for a cosmic-variance-limited galaxy 
survey as in Fig. Q] Here we consider the galaxy survey covering the redshift range 0 < z< z max - The top curve shows the purely geometrical 
constraint on (Eq. [7J. giving the fundamental limit in a sense that the constraint is free of any gravity theory or dark energy model and rests 
on the geometrical nature of the four-dimensional spacetime of the expanding homogeneous and isotropic universe. The cosmic-variance- 
limited survey can achieve a precision better than 0.1% such as cr(Cl K ) - 0.0008. If we assume Einstein gravity (GR) and employ a model of 
dark energy, the curvature constraint can be improved. However, the precision depends on which dark energy model we adopt. If we assume 
the dark energy model parametrized by its equation of state, w DE (z) = w 0 + u> a (l - a) around the fiducial model with the cosmological constant, 
the curvature constraint can be significantly improved. If we restrict ourselves to the cosmological constant (uide = -1), the best curvature 
constraint is achieved (as denoted by the lower curve), but the constraint is model dependent. 


is different from that of the curvature (cT 2 ), the BAO distance 
measurements determine the curvature and the dark energy 
parameters simultaneously. If the dark energy has an (even 
very small) additive contribution, such as dpoE K 1/a 2 , the 
contribution leaves a strong degeneracy with the curvature. 

Nevertheless it would be useful to estimate the accuracy of 
the curvature determination from the cosmic-variance-limited 
BAO measurements when assuming the above dark energy 
model. To do this, we use the method in Ref. l23l to prop¬ 
agate the accuracies of BAO distance measurements into esti¬ 
mations of cosmological parameters. To be more precise, we 
compute the submatrix of the inverted Fisher matrix of CMB 
plus the BAO measurements: 

invert , use submatrix , 

F => [F]- 1 => [F]- m \ n ,. (23) 

The submatrix Fm'N' includes marginalization over other pa¬ 
rameters such as galaxy bias (see Eq. E3- Then the elements 
denoted by M',N' include parameters that determine the dis¬ 


tances among the set of parameters (Eq. 17 1 , and then we 
make the parameter forecasts as follows: 


Pm' — {O m0 , Q m o/r, D A {zd, ^ h (a) } 

_„ _ dp M ’ _ dp M ' 

=> Tde,;; = ~z — r m'N'—z — 

oqi oc[ j 

=> qi = {£2m()/i-, ^DEO, U> 0 , w a ). (24) 

The second line denotes a projection of the Fisher matrix onto 
a new parameter space denoted by q t . Then the third line de¬ 
notes a derivation of the marginalized error on the parameter 
including the curvature and the dark energy parameters. 

The dashed curves in Fig. [2] show the marginalized error of 
the curvature parameter when assuming GR to model the Hub¬ 
ble expansion history and employing the dark energy model, 
around the fiducial cosmological constant model. By restrict¬ 
ing the analysis to a narrower range of dark energy, the curva¬ 
ture constraint can be improved. If assuming the cosmologi¬ 
cal constant model (the ACDM model), the cosmic-variance- 
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FIG. 3: Similar to the previous figure, but this figure shows how the curvature constraints are degraded for a galaxy survey with a finite number 
density of galaxies in each redshift bin. Here we consider the same number density in different redshift bins for simplicity. 


limited accuracy can be as high as cr(Q K ) - 2 x 10 4 . The 
figure shows that the curvature accuracy gets saturated more 
quickly at z ~ 1, compared to the solid curve without any prior 
on dark energy. This is for two reasons. First, if the dark en¬ 
ergy model is employed, a degeneracy between the curvature 
and the dark energy parameters can be well broken by adding 
the BAO measurements in multiple redshifts, where the dark 
energy around the cosmological constraint affects the expan¬ 
sion history only up to z ~ 1. Second, the sound horizon scale 
can be well determined by adding the BAO measurements up 
to z ~ 1 as we stated above. For these reasons, the BAO mea¬ 
surements at z ~ 1 add relatively less information on the cur¬ 
vature when the dark energy models are employed, compared 
to the case without any prior on dark energy (solid curve). 

In reality, the accuracy of the curvature determination is de¬ 
graded by a finite number of sample galaxies for a spectro¬ 
scopic survey. Figure [3] shows how the curvature estimation 
is degraded as a function of the number density of galaxies. 
Here we assume a constant number density of galaxies in each 
redshift bin for simplicity. The planned galaxy surveys (such 
as the Subaru PFS, Euclid and WFIRST) are designed so as 
to satisfy a requirement of nP g (k) i 1 at k ^ 0.1-0.2/z/Mpc, 
corresponding to a number density of h g - 10 5 (/z/Mpc) 3 
Il32l [34 1. In this case, cr(Q K ) - 0.006 or 0.0009 for the pure 
geometry and for the GR+A case, respectively. However, note 


that these future surveys aim at using mainly emission-line 
galaxies such as [Oil] and/or Her emitters as tracers of large- 
scale structure. The emission-line galaxies are a tiny frac¬ 
tion of imaging galaxies at a depth of i ~ 25mag (less than 
1%). By having a wider wavelength coverage (e.g., up to 
infrared wavelengths as proposed by the SPHEREx mission 
B31 ) and/or improving the sensitivity of the spectrograph, a 
high number density of ~ 10 2 (7z/Mpc) 2 for a redshift galaxy 
survey will be feasible in principle. Intergalactic medium sur¬ 
veys, such as Lyman-a forests [e.g., Ref. |46l or radio intensity 
mapping [e.g., Ref.[47] can also be combined with optical and 
infrared redshift surveys. Hence a (nearly) cosmic-variance- 
limited BAO survey with n ~ lCF 2 (/z/Mpc) 3 is not impossible, 
although challenging. 


IV. DISCUSSION 

We have estimated the accuracy of the geometrical deter¬ 
mination of the spatial curvature from the BAO distance mea¬ 
surements. The BAO is unique in the sense that it allows for a 
measurement of the radial distance at the redshift of the galaxy 
survey, while other methods - such as type-la supernova and 
gravitational lensing - measure the luminosity distance which 
is equivalent to the angular diameter distance. Hence, combin- 
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ing the radial and angular diameter distances allows us to con¬ 
strain the spatial curvature without making many assumption 
about the theory of structure formation (other than the BAO 
scale) such as galaxy bias or any model of dark energy. We 
showed that an all-sky, cosmic-variance-limited galaxy survey 
covering up to a high redshift of z ~ 4 allows for a curvature 
determination to an accuracy of cr{Q. K ) - 10 3 . If we as¬ 
sume a simple model of dark energy that is parametrized by 
two equation-of-state parameters, uu{a) — wq + w a (\ - a), it 
allows for an accuracy of cr(D. K ) - a few x 10 1 2 3 4 5 6 7 , although the 
constraint is considered model dependent in this case. 

If there is another method of measuring the Hubble expan¬ 
sion rate at high redshift, it might further enable us to improve 
the curvature constraint based on the method we discussed in 
this paper. For example, if “red-envelope” galaxies - which 
have formed their stellar population at higher redshift at z ~ 2 
- can be used as a cosmic chronometer (as demonstrated in 
Ref. l48l ). the ages of such red galaxies from detailed spec¬ 
troscopic observations can be used to estimate the Hubble ex¬ 
pansion rate at the galaxy redshift. However, this method still 
seems to be limited by astrophysical systematic effects; a fur¬ 
ther careful study needs to be done. The spectroscopic obser¬ 
vations of an extremely large-aperture telescope, in combina¬ 
tion with the optical frequency comb technique, in principle 
allow for a precision measurement of the Hubble expansion 
rate in the high-redshift Universe [e.g.. Ref. [491 ■ 

There is an interesting, physical contamination to the esti¬ 
mation of the curvature. A coherent density contrast across a 
survey region or a local patch in an inhomogenous universe 
can appear as a curved universe even if the global geometry 
is flat (501 ED. As can be found from Eq. (43) in Ref. (52l . 
the induced, apparent curvature is expressed as 0 ! " cal ~ () h 
with an 0(1) prefactor, where c) h is the coherent density con¬ 
trast. If a local universe is embedded into a coherent under¬ 
density region up to z ~ 1, the comoving survey volume is 
estimated as V(< z = 1) ~ 50 (Gpc//t) 3 for an all-sky sur¬ 
vey. The ACDM model predicts <Tb ~ a few x 10 4 for the 
expected rms of the density contrast for the volume, as can be 
found from Fig. 1 in Ref. l53l [also see Ref.EH- Hence such 
a lcrb negative density contrast could degrade a determina¬ 
tion of the global curvature or completely mimic the global 
curvature, if the true curvature is as small as O k ~ 10 4 . 
However, as stressed in Refs. 15314561 such a coherent density 
contrast also causes characteristic, scale-dependent modifica¬ 
tions in the growth of all small-scale stmctures (such as weak 
lensing and cluster abundance) compared to their flat ACDM 


expectations. Hence, although the apparent curvature effects 
themselves are very intriguing to explore, the effects can in 
principle be distinguished by combining the geometrical BAO 
measurements with probes of large-scale structure growth. In 
addition, even if such a coherent underdensity contrast exists, 
it would be very difficult (or there would only be an incredi¬ 
bly small chance) to have a situation that we on the Earth are 
located at the center of such a large-scale void. Hence having 
a wider-area coverage to explore anisotropic curvature effects 
on the sky as well as having a larger redshift leverage to ex¬ 
plore the void boundary should help distinguish the apparent 
curvature effects observationally. 

Constraining the curvature together with properties of the 
primordial scalar perturbations and the primordial gravita¬ 
tional wave is an important direction to explore. If a nonzero 
spatial curvature is found from cosmological observations, it 
will definitely change our view and understanding of the Uni¬ 
verse, so it is worth exploring with any possible means of cur¬ 
rent or future observational data sets. 
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